Use of nucleotide analogs in the analysis of oligonucleotide mixtures and highly multiplexed nucleic acid sequencing

ABSTRACT

Methods kits and systems that use nucleotide analogs to confer increased accuracy and improved resolution in the analysis and sequencing of oligonucleotide mixtures are provided.

RELATED APPLICATIONS

[0001] This application is a divisional of allowed U.S. application Ser. No. 09/880,988, filed Jun. 13, 2001, to Charles R. Cantor and Fouad A. Siddiqi, entitled “USE OF NUCLEOTIDE ANALOGS IN THE ANALYSIS OF OLIGONUCLEOTIDE MIXTURES AND IN HIGHLY MULTIPLEXED NUCLEIC ACID SEQUENCING.” U.S. application Ser. No. 09/880,988, claims the benefit of priority under 35 U.S.C. §119(e) to U.S. provisional application Serial No. 60/211,356, filed Jun. 13, 2000, to Charles R. Cantor and Fouad A. Siddiqi, entitled USE OF NUCLEOTIDE ANALOGS IN THE ANALYSIS OF OLIGONUCLEOTIDE MIXTURES AND IN HIGHLY MULTIPLEXED NUCLEIC ACID SEQUENCING.” Accordingly, the benefit of the priority date of U.S. provisional application Serial No. 60/211,356 and U.S. application Ser. No. 09/880,988 is claimed. The subject matter of each of U.S. provisional application Serial No. 60/211,356 and U.S. application Ser. No. 09/880,988 is incorporated by reference in it entirety.

[0002] Subject matter described herein was developed under NSF Grant No. Ger-9452651. The Government can have certain rights therein.

FIELD OF THE INVENTION

[0003] This invention relates to methods, particularly mass spectrometric methods, for the analysis and sequencing of nucleic acid molecules.

DESCRIPTION OF THE BACKGROUND

[0004] Since the recognition of nucleic acid as the carrier of the genetic code, a great deal of interest has centered around determining the sequence of that code in the many forms in which it occurs. Two studies made the process of nucleic acid sequencing, at least with DNA, a common and relatively rapid procedure practiced in most laboratories. The first describes a process whereby terminally labeled DNA molecules are chemically cleaved at single base repetitions (A. M. Maxam and W. Gilbert, Proc. Natl. Acad. Sci. USA 74:560-64, 1977). Each base position in the nucleic acid sequence is then determined from the molecular weights of fragments produced by partial cleavage. Individual reactions were devised to cleave preferentially at guanine, at adenine, at cytosine and thymine, and at cytosine alone. When the products of these four reactions are resolved by molecular weight, using, for example, polyacrylamide gel electrophoresis, DNA sequences can be read from the pattern of fragments on the resolved gel.

[0005] In another method DNA is sequenced using a variation of the plus-minus method (Sanger et al. (1977) Proc. Natl. Acad. Sci. USA 74:5463-67, 1977). This procedure takes advantage of the chain terminating ability of dideoxynucleoside triphosphates (ddNTPs) and the ability of DNA polymerase to incorporate ddNTPs with nearly equal fidelity as the natural substrate of DNA polymerase, deoxynucleoside triphosphates (dNTPs). Briefly, a primer, usually an oligonucleotide, and a template DNA are incubated in the presence of a useful concentration of all four dNTPs plus a limited amount of a single ddNTP. The DNA polymerase occasionally incorporates a dideoxynucleotide that terminates chain extension. Because the dideoxynucleotide has no 3′-hydroxyl, the initiation point for the polymerase enzyme is lost. Polymerization produces a mixture of fragments of varied sizes, all having identical 3′ termini. Fractionation of the mixture by, for example, polyacrylamide gel electrophoresis, produces a pattern that indicates the presence and position of each base in the nucleic acid. Reactions with each of the four ddNTPs permits the nucleic acid sequence to be read from a resolved gel.

[0006] These procedures are cumbersome and are limited to sequencing DNA. In addition, with conventional procedures, individual sequences are separated by, for example, electrophoresis using capillary or slab gels, which slow. Mass spectrometry has been adapted and used for sequencing and detection of nucleic acid molecules (see, e.g., U.S. Pat. Nos. (6,194,144; 6,225,450; 5,691,141; 5,547,835; 6,238,871; 5,605,798; 6,043,031; 6,197,498; 6,235,478; 6,221,601; 6,221,605). In particular, Matrix-Assisted Laser Desorption/lonization (MALDI) and ElectroSpray Ionization (ESI), which allow intact ionization, detection and exact mass determination of large molecules, i.e. well exceeding 300 kDa in mass have been used for sequencing of nucleic acid molecules.

[0007] A further refinement in mass spectrometric analysis of high molecular weight molecules was the development of time of flight mass spectrometry (TOF-MS) with matrix-assisted laser desorption ionization (MALDI). This process involves placing the sample into a matrix that contains molecules that assist in the desorption process by absorbing energy at the frequency used to desorb the sample. Time of flight analysis uses the travel time or flight time of the various ionic species as an accurate indicator of molecular mass. Due to its speed and high resolution, time-of-flight mass spectrometry is well-suited to the task of short-range, i.e., less than 30 base sequencing of nucleic acids. Since each of the four naturally occurring nucleotide bases dC, dT, dA and dG, also referred to herein as C, T, A and G, in DNA has a different molecular weight,

M_(c)=289.2

M_(T)=304.2

M_(A)=313.2

M_(G)=329.2,

[0008] where M_(c), M_(T), M_(A), M_(G) are average molecular weights in daltons of the nucleotide bases deoxycytidine, thymidine, deoxyadenosine, and deoxyguanosine, respectively, it is possible to read an entire sequence in a single mass spectrum. If a single spectrum is used to analyze the products of a conventional Sanger sequencing reaction, where chain termination is achieved at every base position by the incorporation of dideoxynucleotides, a base sequence can be determined by calculation of the mass differences between adjacent peaks. In addition, the method can be used to determine the masses, lengths and base compositions of mixtures of oligonucleotides and to detect target oligonucleotides based upon molecular weight.

[0009] MALDI-TOF mass spectrometry for sequencing DNA using mass modification (see, e.g., U.S. Pat. Nos. 5,547,835, 6,194,144; 6,225,450; 5,691,141 and 6,238,871) to increase mass resolution is available. The methods employ conventional Sanger sequencing reactions with each of the four dideoxynucleotides. In addition, for example for multiplexing, two of the four natural bases are replaced; dG is substituted with 7-deaza-dG and dA with 7-deaza-dA.

[0010] A variety of techniques and combinations thereof have been directed to improving the level of accuracy in determining the nucleotide compositions of mixtures of oligonucleotides using mass spectrometry, and many of these methods employ nucleotide analogs. For example, Muddiman et al. (Anal. Chem., 69(8): 1543-1549, 1997) discusses an algorithm for the unique definition of the base composition of PCR-amplified products, especially longer (>100 bp) oligonucleotides. The algorithm places a constraint on the otherwise large number of possible base compositions for long oligonucleotides by taking into account only those masses (measured by electrospray ionization mass spectrometry) that are consistent with that of their denatured complementary strands, assuming Watson-Crick base-pairing. In addition, the algorithm imposes the constraint of known primer compositions, since the primer sequences are known, and this constraint becomes especially significant with shorter PCR products whose mass of “unknown” sequence relative to that of the primer mass is small. Muddiman et al. also discusses invoking additional measurements for defining the base composition with even greater accuracy. These include the possibility of post-modifying the PCR product using e.g., dimethyl sulfate to selectively methylate every “G” in the PCR product, or using a modified base during PCR amplification, conducting mass measurements on the modified oligonucleotides, and comparing the mass measurements with those of the unmodified complementary strands.

[0011] Chen et al. (Anal. Chem., 71 (15): 3118-3125, 1999) reports a method that combines stable isotope ¹³C/¹⁵N labelling of PCR products with analysis of the mass shifts by MALDI-TOF mass spectrometry. The mass shift due to labelling of a single type of nucleotide (i.e, A, T, G or C) reveals the number of that type of nucleotide in a given fragment. While the method is useful in the measurement and comparison of nucleotide compositions of homologous sequences for sequence validation and in scoring polymorphisms, tedious repetitive sequencing reactions (using the four different labelled nucleotides) and mass spectrometric measurements are required.

[0012] Hence there is a need in the art for methods that (i) unambiguously assign nucleotides in a sequence, and, (ii) resolve large numbers of oligonucleotides that have the same length, different base compositions, and nearly equal (i.e., less than or equal to about 1 dalton difference) molecular weights. Therefore it is an object herein to provide methods that solve such problems

SUMMARY OF THE INVENTION

[0013] Provided herein are methods for sequencing and detecting nucleic acids using techniques, such as mass spectrometry and gel electrophoresis, that are based upon molecular mass. The methods use deoxynucleotide analogs, modified nucleotide terminators and/or mass-labeled primers in one or more reactions for sequencing or detection protocols that involve primer extension, and analyze these products from more than one oligonucleotide on, for example, a single mass spectrum. This provides a means for accurate detection and/or sequencing of a an oligonucleotide and is particularly advantageous for detecting or sequencing a plurality target nucleic acid molecules in a single reaction using any technique that distinguishes products based upon molecular weight. The methods herein are particularly adapted for mass spectrometric analyses.

[0014] For example, a sequencing method provided herein uses deoxynucleotide analogs, modified nucleotide terminators and/or mass-labeled primers in one or more Sanger sequencing reactions, and analyzes these products from more than one oligonucleotide on a single mass spectrum. In particular, a plurality of primers can be used to simultaneously sequence a plurality of nucleic acid molecules or portions of the same molecule. By extending the primers with mass-matched nucleotides, the resulting products mass shifts that are periodically related to the size of the original primer.

[0015] As a result, the sequence of any given oligonucleotide can be determined with a high level of accuracy, and also mixtures of a number of sequences can be multiplexed in a single mass spectrum. The limit on the number of oligonucleotides that can be sequenced simultaneously is governed by the base periodicity, the maximum mass shift, and the resolving power of analytical tool, such as the mass spectrometer. The base periodicity and maximum mass shift can be carefully engineered for optimal resolution and accuracy, depending on the number of sequences to be simultaneously analyzed, and the information desired; as many sequences as desired can be sequenced simultaneously especially in the detection and scoring of single nucleotide polymorphisms, insertions, deletions and other mutations.

[0016] In another embodiment, a target nucleic acid molecule is sequenced using mass-matched nucleotides and chain terminating nucleotides. For example, a primer is annealed to a target nucleic acid, the primer is extended in the presence of chain-terminating nucleotides and mass-matched nucleotides to produce extension products, the masses of the extension products follow a periodic distribution that is determined by the mass of the mass-matched nucleotides, and the sequence of the target nucleic acid is determined from the mass shift of each extension product from its corresponding periodic reference mass by virtue of incorporation of the chain terminator. The mass-matched nucleotides all have identical masses, and each chain terminating nucleotide has a distinct mass that differs from that of the other chain terminating nucleotides. This results in unique predetermined values of mass shift corresponding to each chain terminating nucleotide and based upon the original primer.

[0017] This method is adaptable for any sequencing method or detection method that relies upon or includes chain extension. These methods include, but are not limited to, sequencing methods based upon Sanger sequencing, and detection methods, such as primer oligo base extension (PROBE) (see, e.g., U.S. application Ser. No. 6,043,031; allowed U.S. application Ser. No. 09/287,679; and 6,235,478), that rely include a step of chain extension.

[0018] Also, contemplated are methods, such as haplotyping methods, in which two mutations in the same gene are detection are provided. A detector (primer) oligonucleotide is to the hybridized to the first mutation and the primer is extended with mass-matched nucleotides and appropriately selected chain terminator(s) to detect the second mutation.

[0019] In other embodiments, a plurality of target nucleic acids can be multiplexed in a single reaction measurement by annealing each target nucleic acid to a primer of distinct molecular weight each primer is then extended with mass-matched nucleotides and chain terminators in formats that depend upon whether detection or sequencing is desired. These methods are particularly useful for methods of detection in which a primer is hybridized to a plurality of target nucleic acid molecules, such as immobilized nucleic acid molecules, hybrids separated from unhybridized nucleic acids and the detectors detected. Such methods include PROBE, in which case the extension reaction is performed in the presence chain terminators and mass matched deoxynucleotides.

[0020] The primers of distinct molecular weight can be selected to differ in molecular weight by a value that is greater than the maximum mass shift, i.e., the difference in molecular weight between the heaviest and the lightest nucleotide terminators in chain extension reactions. The difference in molecular weight between the primers for a plurality of target nucleic acids can be selected to be least 20 daltons greater than the maximum mass shift to account for the finite band width of the peaks.

[0021] The number of molecules that can be multiplexed is governed by the periodicity, the maximum mass shift, and the resolving power of the sequence detection instrument. In some embodiments, about 7 to about 25 or more molecules can be multiplexed. For scoring single nucleotide polymorphisms, only a single nucleotide terminator is required (depending on the base identity of the single nucleotide polymorphism). In this case, the maximum mass shift required is identically zero, so that larger numbers of molecules, greater than 25, 35, 50 and more, can be multiplexed, depending on the resolving power of the sequencing format, and for mass spectrometry the instrument. Depending on the amount of sequence information desired, one, two or three rather than four types of nucleotide terminators (corresponding to each of the four nucleic acid bases) can be used.

[0022] In other embodiments, the mass shift is obtained using pair-matched nucleotides, i.e., the mass of each nucleotide base-pair is selected so that the masses of all pairs are identical. In one embodiment thereof, the following steps are performed: (i) the target nucleic acid is copied or amplified by a method such as PCR in the presence of the pair-matched nucleotide set prior to the sequencing or detection reaction; (ii) the target nucleic acid is denatured, and a partially duplex hairpin primer is annealed and ligated to the single-stranded template; (iii) the primer is extended in the presence of chain terminating nucleotides and pair-matched nucleotides to produce extension products, where the masses of the extension products follow a periodic distribution that is determined by the mass of the pair-matched nucleotide set, and, (iv) the target nucleic acid is detected by virtue of its molecular weight or its sequence is determined from the mass shift of each extension product from its corresponding periodic reference mass.

[0023] In another embodiment, the mass of each terminating base pair is unique and resolvable, so that the mass shifts corresponding to each terminating base pair are unique. The nucleotide terminators are optionally mass-matched or can be of distinct masses as long as distinct values of mass shift are obtained for each terminating base pair.

[0024] In another embodiment, the extension products are treated to produce blunt-ended double-stranded extension products by methods known to those of skill in the art, such as the use of single-strand specific nucleases. In an aspect of this embodiment, a plurality of target nucleic acids can be multiplexed in a single reaction by annealing each target nucleic acid to a primer of distinct molecular weight. The primers can be selected to differ in molecular weight by a value that is greater than the maximum mass shift, i.e., the difference in molecular weight between the heaviest and the lightest nucleotide terminating base pairs. Since double stranded nucleic acid can be analyzed, the effective sequence read is halved relative to the embodiment employing mass-matched nucleotides, but the number of molecules that can be multiplexed is doubled, due to the increase in period (the value of the mass of a base pair, rather than a single mass-matched nucleotide). In exemplary embodiments, about 14 to about 50 sequences are multiplexed. In detection embodiments, about 50 or more molecules can be simultaneously multiplexed since only a single terminating base pair is added in the extension reaction.

[0025] In another embodiment, the chain termination reactions are carried out separately using a standard nucleotide terminator, pair-matched nucleotides, and mass-labeled primers, if modified nucleotide terminators which are either mass-matched or provide distinct values of mass shift for each terminating base pair are not available. The reactions are pooled prior to detection or sequence analysis. In one embodiment, the mass-labeled primers can have distinct values of molecular weight that give rise to unique values of mass shift or positional mass difference for each terminating base.

[0026] In and other method provided herein, a population of nucleic acids having the same length but different base compositions can be resolved by synthesizing the nucleic acids in the presence of a nucleotide analog to produce synthesized nucleic acids having incorporated the nucleotide analog, where the nucleotide analog is selected to optimally separate the masses of the population of nucleic acids according to their individual base compositions. For example, the nucleotide analog or analogs are selected to separate the population of nucleic acids according to base composition by greater than 1 dalton. In another embodiment, the nucleotide analog or analogs are selected to separate the population of nucleic acids according to base composition by mass values of about 3 daltons to about 8 daltons, depending on the choice of analog and on the resolving power of the detection instrument. In other embodiments, the nucleotide analog or analogs can be selected to restrict oligonucleotides having the same length to have the same mass, i.e., a peak separation of zero, regardless of differences in base composition, such as in detection methods, where it is desirable to separate populations of oligonucleotides according to their length.

[0027] Nucleic acid molecules that contain mass-matched nucleotides and/or pair-matched nucleotides are provided.

[0028] Also provided are combinations for practicing the methods provided herein. For instance, in one embodiment, the combinations include a set of mass-matched deoxynucleotides. In another embodiment, the combinations a set of pair-matched nucleotides and a set of mass-matched chain terminating nucleotides. In another embodiment, the combination includes a set of pair-matched nucleotides and chain terminating nucleotides which form terminating base pairs of distinctly different molecular weight. In yet another embodiment, the combination includes a set of pair-matched nucleotides and mass-labeled primers. In other embodiments, mass-staggered primers can be added to as optional components.

[0029] Kits containing the combinations with optional instructions and/or additional reagents are also provided. The kits contain the reagents as described herein and optionally any other reagents required to perform the reactions. Such reagents and compositions are packaged in standard packaging known to those of skill in the art. Additional vials, containers, pipets, syringes and other products for sequencing can also be included. Instructions for performing the reactions can be included.

[0030] Also provided herein are methods for optimization of the analysis of base compositions of mixtures of oligonucleotides by mass spectrometry. A single spectrum can be used to resolve a very large number of oligonucleotides having the same length but different molecular weights by incorporating a nucleotide analog into the oligonucleotides in the mixture such that the peaks are no closer than a minimum value called peak separation. The peak separation can be tailored by careful selection of the nucleotide analog and of a mass spectrometer with the desired resolving power.

[0031] The methods herein permit unambiguous and accurate analysis of the sequences or molecular weights of large numbers of oligonucleotides in a single mass spectrum by combining the rapidity of mass spectrometry with the resolving power of nucleotide analogs which are carefully selected and incorporated into the oligonucleotide mixture according to the desired application.

[0032] Other features and advantages will be apparent from the following detailed description and claims.

BRIEF DESCRIPTION OF THE FIGURES

[0033]FIG. 1 shows that when a single spectrum is used to analyze the products of a conventional Sanger sequencing reaction, where chain termination is achieved at every base position by the incorporation of dideoxynucleotides, the base sequence can be determined by calculation of the mass differences between adjacent peaks (FIGS. 1a and 1 b).

[0034]FIG. 2 shows implementation of forced mass modulation using mass-matched deoxynucleotides. FIG. 2a is a simulated mass spectrum showing the products and molecular masses of a reaction carried out with a suitable polymerase in the presence of a mass-matched nucleotide set (“dN”) and the four standard dideoxynucleotide terminators. The base periodicity is the mass of dN, or 310 daltons. FIG. 2b shows a target second sequence resolved on the same mass spectrum shown in FIG. 2a, using a primer heavier by 77 daltons. The peaks corresponding to the reaction products from the first target sequence can fall within the spectrum in FIG. 2b, which can never intersect peaks from the second target sequence. This permits unambiguous resolution of both sequences each peak can be uniquely assigned to a nucleotide, a base position, and a target sequence.

[0035]FIG. 3 shows four different sequences resolved in a single spectrum using a set of mass-staggered primers that are separated in mass by integer multiples of 77 daltons (77, 154, and 231 daltons).

[0036]FIG. 4 shows the general implementation of a forced mass modulation method using pair-matched nucleotides, for the analysis of sequencing reaction products as double-stranded structures. The steps in the reaction are as follows: a) a partially duplex hairpin primer with a 3′ overhang and a 5′ phosphate group is annealed and ligated to the single stranded target sequence; b) the resulting partially duplex structure is subjected to a sequencing reaction using the pair-matched nucleotide set described above along with the set of mass-matched terminators (ddM); c) products resulting from sequencing reaction b); and, d) the products c) from the sequencing reaction are exposed to a strict single strand-specific nuclease that results in the production of blunt-ended hairpin structures ready for analysis by mass spectrometry.

[0037]FIG. 5 shows the products and molecular masses of the nuclease digestion elucidated in FIG. 4d, along with a simulated mass spectrum.

[0038]FIG. 6 shows three sequence variants (FIG. 6a) that differ from each other only at a single base position sequenced by a conventional Sanger reaction. FIG. 6b is a simulated mass spectrum of all reaction products shown in FIG. 6a. FIG. 6c is a graph representing the valid sequence permutations that can be elucidated from the mass spectrum shown in FIG. 6b. Boxed values are fragment masses, solid arrows show valid sequence branches, dashed arrows represent spurious branches. In practice valid branches are indistinguishable from spurious ones. FIG. 6d is a set of sequences consistent with the graph shown in FIG. 6c. Spurious sequence reconstructions are shown in lowercase letters, valid ones in uppercase letters.

[0039]FIG. 7a shows the three related sequence variants from FIG. 6 sequenced by Forced Mass Modulation using a single primer and the mass-matched nucleotide set from FIG. 2 with the standard dideoxy terminators. The positions of the differing bases are shown by solid arrows. Reaction products are shown along with their respective molecular masses. Reaction products of variant #2 whose masses differ from those of variant #1 are marked with by (*). Reaction products of variant #3 whose masses differ from those of variant #1 are marked by (**). FIG. 7b is a simulated mass spectrum of all reaction products shown in FIG. 7a along with a sequence graph. The shaded regions represent the only valid mass ranges that can be assumed by the reaction products from FIG. 7a. The base periodicity is 310 daltons. FIG. 7c is a consensus sequence derived from the data shown in FIG. 7b. FIG. 7d is an expansion of the consensus sequence shown in FIG. 7c. Spurious reconstructions are shown in lowercase letters, valid ones in uppercase letters. Note that there is only a single spurious reconstruction, as opposed to the eleven errant sequences reconstructed from the Sanger reaction described in FIG. 6.

[0040]FIG. 8 shows the base composition density distributions for the total set of possible 7-base oligonucleotides using three different nucleotide sets. Note that for the set of naturally occurring bases, nearly every base composition has its own distinct mass value, but most of these mass values are spaced only one dalton from each other. Increasing the peak separation also markedly increases the average number of base compositions per observed mass, particularly for those masses in the center of the range.

DETAILED DESCRIPTION OF THE INVENTION

[0041] Definitions

[0042] Unless defined otherwise, all technical and scientific terms used herein have the same meaning as is commonly understood by one of skill in the art to which this invention belongs. All patents, patent applications, Genbank and other sequence repository sequences, and publications referred to herein are incorporated by reference.

[0043] As used herein, a biopolymer includes, but is not limited to, nucleic acid, proteins, polysaccharides, lipids and other macromolecules. Nucleic acids include DNA, RNA, and fragments thereof. Nucleic acids may be derived from genomic DNA, RNA, mitochondrial nucleic acid, chloroplast nucleic acid and other organelles with separate genetic material.

[0044] As used herein “nucleic acid” refers to polynucleotides such as deoxyribonucleic acid (DNA) and ribonucleic acid (RNA). The term should also be understood to include, as equivalents, derivatives, variants and analogs of either RNA or DNA made from nucleotide analogs, single (sense or antisense) and double-stranded polynucleotides. Deoxyribonucleotides include deoxyadenosine, deoxycytidine, deoxyguanosine and deoxythymidine. For RNA, the uracil base is uridine.

[0045] As used herein, “forced mass modulation” refers to methods provided herein that use deoxynucleotide analogs, modified nucleotide terminators, mass-labeled primers, mass-staggered primers and other such nucleotides, nucleic acids and analogs thereof, to unambiguously assign peak positions of mass fragments of oligonucleotides according to their base position, base identity, and target sequence from which the fragments arose. The method is used to sequence, detect or identify single oligonucleotide or plurality thereof. Hence the method is used, for example for muliplex sequencing and detection of nucleic acid molecules among mixtures thereof.

[0046] As used herein, “nucleotides” include, but are not limited to, the naturally occurring nucleoside mono-, di-, and triphosphates: deoxyadenosine mono-, di- and triphosphate; deoxyguanosine mono-, di- and triphosphate; deoxythymidine mono-, di- and triphosphate; and deoxycytidine mono-, di- and triphosphate (referred to herein as dA, dG, dT and dC or A, G, T and C, respectively). Nucleotides also include, but are not limited to, modified nucleotides and nucleotide analogs such as deazapurine nucleotides, e.g., 7-deaza-deoxyguanosine (7-deaza-dG) and 7-deaza-deoxyadenosine (7-deaza-dA) mono-, di- and triphosphates, deutero-deoxythymidine (deutero-dT) mon-, di- and triphosphates, methylated nucleotides e.g., 5-methyldeoxycytidine triphosphate, ¹³C/¹⁵N labelled nucleotides and deoxyinosine mono-, di- and triphosphate. For those skilled in the art, it will be clear that modified nucleotides and nucleotide analogs can be obtained using a variety of combinations of functionality and attachment positions.

[0047] As used herein, a complete set of chain-elongating nucleotides refers to the four different nucleotides or analogs thereof that hybridize to each of the four different bases comprising the nucleic acid template.

[0048] As used herein, the term “mass-matched nucleotides” refers to a set of nucleotide analogs wherein each analog is of identical mass to each of the other analogs. For example, analogs of dA, dG, dC and dT can form a mass-matched nucleotide set, when each analog is selected to have the same molecular weight as the others in the set. Mass-matched nucleotide sets can be identified by selecting chemically modified derivatives of natural bases or by the use of a universal base analog such as deoxyinosine or 5-nitroindole and 3-nitropyrrole (5-nitroindole and 3-intropyrrole can be in the dideoxy form) which can form base pairs with more than one of the natural bases. Others include, 3-methyl 7-propynyl isocarbostyril, 5-methyl isocarbostyril, and 3-methyl isocarbostyril. As a result, oligonucleotides that contain such bases differ in molecular weight only as a function of length thereof. Furthermore, incorporation of a single nucleotide(s) that is (are) not in the set renders such the oligonucleotide(s) readily identifiable by mass, particularly by spectrometric analysis.

[0049] As used herein, the term “pair-matched nucleotides” refers to a nucleotide set in which the nucleotide analogs are selected such that the total mass each base pair is identical. For example, replacing dG with the nucleotide analog 7-deaza-dG forces the mass of each base pair, i.e., (dA+dT) and (dC+7-deaza-dG) to be identical. Exemplary pair-matched nucleotides, include, but are not limited to, 7-deaza-dA+phosphorothioate-dT ((312.2+320.2)=632.4 Da) and 5-methyl-dC+dG ((303.2+329.2)=632.4 Da); phosphorothioate-7-deaza-dA+dU ((328.2+290.2)=618.4 Da) and dC+dG=((289.2+329.2)=618.4 Da), and other such pairs that may be readily selected. Another exemplary set of mass-matched nucleotides with a molecular mass of 328.2: 7-deaza-dG, phosphorothioate-7-deaza-dA, 5-propynyl-dU and 5-cyanomethyl-2′-deoxycytidine.

[0050] As used herein, the term “nucleotide terminator” or “chain terminating nucleotide” refers to a nucleotide analog that terminates nucleic acid polymer (chain) extension during procedures wherein a DNA template is being sequenced or replicated. The standard chain terminating nucleotides, i.e., nucleotide terminators include 2′,3′-dideoxynucleotides (ddATP, ddGTP, ddCTP and ddTTP, also referred to herein as dideoxynucleotide terminators). As used herein, dideoxynucleotide terminators also include analogs of the standard dideoxynucleotide terminators, e.g., 5-bromo-dideoxyuridine, 5-methyl-dideoxycytidine and dideoxyinosine are analogs of ddTTP, ddCTP and ddGTP, respectively.

[0051] As used herein, “mass-matched terminators” refers to a set of nucleotide terminators that are selected such that each analog of ddA, ddG, ddC and ddT making up the mass-matched set has exactly the same molecular weight. Mass-matched terminator sets can be constructed by selecting chemically modified derivatives of standard dideoxynucleotides or by the use of a universal dideoxynucleotide analog that form base pairs with more than one of the natural bases. Exemplary mass-matched nucleotides include, but are not limited to, 3-methyl 7-propynyl isocarbostyril, 5-methyl isocarbostyril and 3-methyl isocarbostyril. As used herein, the terms “oligonucleotide” or “nucleic acid” refer to single-stranded and/or double-stranded polynucleotides such as deoxyribonucleic acid (DNA), and ribonucleic acid (RNA) as well as derivatives of either RNA or DNA into which nucleotide or dideoxynucleotide analogs have been incorporated. Also included in the term “nucleic acid” are analogs of nucleic acids such as peptide nucleic acid (PNA), phosphorothioate DNA, and other such analogs and derivatives.

[0052] As used herein, “nucleotide composition” or “base composition” refers to the numerical ratio of the four nucleotide bases relative to each other in an oligonucleotide.

[0053] As used herein, a target nucleic acid refers to any nucleic acid of interest in a sample. It can contain one or more nucleotides. A target nucleotide sequence refers to a particular sequence of nucleotides in a target nucleic acid molecule. Detection or identification of such sequence results in detection of the target and can indicate the presence or absence of a particular mutation or polymorphism.

[0054] As used herein, “partially duplex hairpin” refers to a partially self-complementary oligonucleotide, which forms intramolecular base-pairs within its self-complementary region, leaving a “loop” of bases at one end of the molecule and a single-stranded “overhang” region at the other end. Thus, the oligonucleotide assumes a hairpin-like motif. “Blunt-ended hairpin structures”, as referred to herein, are similar to the partially duplex hairpin structures with the exception that they do not have a single-stranded “overhang” region.

[0055] As used herein, “base periodicity” or “period” (P_(base)) refers to the quasi-periodic distribution of the molecular weights of products obtained using Forced Mass Modulation. The base periodicity results from either the mass of the mass-matched deoxynucleotide set, or the mass of the pair-matched deoxynucleotide set, or from the modified chain terminators depending on the embodiment implemented. The base sequence or nucleic acid molecule identity is encoded in the pattern (or detectable therein) in which the observed mass distribution deviates from absolute regular periodicity.

[0056] As used herein, the “periodic reference mass” at base position “n” in any given oligonucleotide molecule, M_(PR)[n], is defined as the sum of: (i) the mass of the primer (M_(primer)) used to sequence the DNA template using Forced Mass Modulation, (ii) the mass of the lightest nucleotide terminator (M_(light)), and, (iii) (n−1) multiple of the base periodicity P_(base).

[0057] As used herein, the “positional mass difference” or “mass shift” at base position “n” in any given oligonucleotide molecule, M_(diff)[n], is defined as the distance in daltons between the observed peak, M_(obs)[n], and the nth periodic reference mass.

[0058] As used herein, the “maximum mass shift” S_(max) is the maximum possible value of the positional mass difference, depending on the choice of mass-matched nucleotides and nucleotide terminators used in the implementation of Forced Mass Modulation. Accordingly, the maximum mass shift can be modulated by the choice of mass-matched nucleotides and nucleotide terminators.

[0059] As used herein, a “primer” refers to an oligonucleotide that is suitable for hybridizing, chain extension, amplification and sequencing. Similarly, a probe is a primer used for hybridization. The primer refers to a nucleic acid that is of low enough mass, typically about between about 5 and 200 nucleotides, generally about 70 nucleotides or less than 70, and of sufficient size to be conveniently used in the methods of amplification and methods of detection and sequencing provided herein. These primers include, but are not limited to, primers for detection and sequencing of nucleic acids, which require a sufficient number nucleotides to form a stable duplex, typically about 6-30 nucleotides, about 10-25 nucleotides and/or about 12-20 nucleotides. Thus, for purposes herein, a primer is a sequence of nucleotides of any suitable length, typically containing about 6-70 nucleotides, 12-70 nucleotides or greater than about 14 to an upper limit of about 70 nucleotides, depending upon sequence and application of the primer.

[0060] As used herein, the term “mass-labeled primers” refers to a set of primers that differ in mass by values that provide distinct and resolvable positional mass differences for each of the four termination reactions in an embodiment of Forced Mass Modulation. In this particular embodiment of Forced Mass Modulation, each of the termination reactions for a given oligonucleotide is carried out separately using each of the mass-labeled primers, and the reaction products are combined prior to obtaining a mass spectrum.

[0061] As used herein, the term “mass-staggered primers” refers to the mass difference (“staggering” of the masses) between the primers used in multiplexed sequencing using Forced Mass Modulation. For resolution of multiple sequences using this method, the differences between the masses of the primers should at least be equal to the maximum mass shift, and is generally greater than the maximum mass shift by at least 20 daltons to account for the finite width of each observed peak.

[0062] As used herein, reference to mass spectrometry encompasses any suitable mass spectrometric format known to those of skill in the art. Such formats include, but are not limited to, Matrix-Assisted Laser Desorption/lonization, Time-of-Flight (MALDI-TOF), Electrospray (ES), IR-MALDI (see, e.g., published International PCT application No. 99/57318 and U.S. Pat. No. 5,118,937), Ion Cyclotron Resonance (ICR), Fourier Transform and combinations thereof. MALDI, particular UV and IR, are among the preferred formats.

[0063] As used herein, mass spectrum refers to the presentation of data obtained from analyzing a biopolymer or fragment thereof by mass spectrometry either graphically or encoded numerically.

[0064] As used herein, pattern with reference to a mass spectrum or mass spectrometric analyses, refers to a characteristic distribution and number of signals (such peaks or digital representations thereof).

[0065] As used herein, signal in the context of a mass spectrum and analysis thereof refers to the output data, which the number or relative number of molecules having a particular mass. Signals include “peaks” and digital representations thereof.

[0066] As used herein, “mass spectrum division multiplexing” is an embodiment of Forced Mass Modulation in which unambiguous resolution of multiple sequences in a single spectrum is possible by judicious selection of mass staggered primers.

[0067] As used herein, “analysis” refers to the determination of certain properties of a single oligonucleotide, or of mixtures of oligonucleotides. These properties include, but are not limited to, the nucleotide composition and complete sequence of an oligonucleotide or of mixtures of oligonucleotides, the existence of single nucleotide polymorphisms between more than one oligonucleotide, the masses and the lengths of oligonucleotides and the presence of a molecule or sequence within molecule in a sample.

[0068] As used herein, “multiplexing” refers to the simultaneous determination of more than one oligonucleotide molecule, or the simultaneous analysis of more than one oligonucleotide, in a single mass spectrometric or other sequence measurement, i.e., a single mass spectrum or other method of reading sequence.

[0069] As used herein, “polymorphisms” refer to variants of a gene or an oligonucleotide molecule that differ at more than one base position. In “single nucleotide polymorphisms”, the variants differ at only a single base position.

[0070] As used herein, amplifying refers to means for increasing the amount of a bipolymer, especially nucleic acids. Based on the 5′ and 3′ primers that are chosen, amplification also serves to restrict and define the region of the genome which is subject to analysis. Amplification can be by any means known to those skilled in the art, including use of the polymerase chain reaction (PCR) etc. Amplification, e.g., PCR must be done quantitatively when the frequency of polymorphism is required to be determined.

[0071] As used herein, “polymorphism” refers to the coexistence of more than one form of a gene or portion thereof. A portion of a gene of which there are at least two different forms, i.e., two different nucleotide sequences, is referred to as a “polymorphic region of a gene”. A polymorphic region can be a single nucleotide, the identity of which differs in different alleles. A polymorphic region can also be several nucleotides in length. Thus, a polymorphism, e.g. genetic variation, refers to a variation in the sequence of a gene in the genome amongst a population, such as allelic variations and other variations that arise or are observed. Thus, a polymorphism refers to the occurrence of two or more genetically determined alternative sequences or alleles in a population. These differences can occur in coding and non-coding portions of the genome, and can be manifested or detected as differences in nucleic acid sequences, gene expression, including, for example transcription, processing, translation, transport, protein processing, trafficking, DNA synthesis, expressed proteins, other gene products or products of biochemical pathways or in post-translational modifications and any other differences manifested amongst members of a population. A single nucleotide polymorphism (SNP) refers to a polymorphism that arises as the result of a single base change, such as an insertion, deletion or change in a base.

[0072] A polymorphic marker or site is the locus at which divergence occurs. Such site may be as small as one base pair (an SNP). Polymorphic markers include, but are not limited to, restriction fragment length polymorphisms, variable number of tandem repeats (VNTR's), hypervariable regions, minisatellites, dinucleotide repeats, trinucleotide repeats, tetranucleotide repeats and other repeating patterns, simple sequence repeats and insertional elements, such as Alu. Polymorphic forms also are manifested as different mendelian alleles for a gene. Polymorphisms may be observed by differences in proteins, protein modifications, RNA expression modification, DNA and RNA methylation, regulatory factors that alter gene expression and DNA replication, and any other manifestation of alterations in genomic nucleic acid or organelle nucleic acids.

[0073] As used herein, “polymorphic gene” refers to a gene having at least one polymorphic region.

[0074] As used herein, “allele”, which is used interchangeably herein with “allelic variant” refers to alternative forms of a gene or portions thereof. Alleles occupy the same locus or position on homologous chromosomes. When a subject has two identical alleles of a gene, the subject is said to be homozygous for the gene or allele. When a subject has two different alleles of a gene, the subject is said to be heterozygous for the gene. Alleles of a specific gene can differ from each other in a single nucleotide, or several nucleotides, and can include substitutions, deletions, and insertions of nucleotides. An allele of a gene can also be a form of a gene containing a mutation.

[0075] As used herein, “predominant allele” refers to an allele that is represented in the greatest frequency for a given population. The allele or alleles that are present in lesser frequency are referred to as allelic variants.

[0076] As used herein, a subject, includes, but is not limited to, animals, plants, bacteria, viruses, parasites and any other organism or entity that has nucleic acid. Among subjects are mammals, preferably, although not necessarily, humans. A patient refers to a subject afflicted with a disease or disorder.

[0077] As used herein, a phenotype refers to a set of parameters that includes any distinguishable trait of an organism. A phenotype can be physical traits and can be, in instances in which the subject is an animal, a mental trait, such as emotional traits.

[0078] As used herein, “resolving power” of a mass spectrometer is the ion separation power of the instrument, i.e., it is a measure of the ability of the mass spectrometer to separate peaks representing different masses. The resolving power R is defined as m/Δm, where m is the ion mass and Δm is the difference in mass between two resolvable peaks in a mass spectrum.

[0079] As used herein, “assignment”¹ refers to a determination that the position of a nucleic acid fragment indicates a particular molecular weight and a particular terminal nucleotide.

[0080] As used herein, “a” refers to one or more.

[0081] As used herein, “plurality” refers to two or more, up to an amount that is governed by the base periodicity, the maximum mass shift, and the resolving power of the mass spectrometer.

[0082] As used herein, an array refers to a pattern produced by three or more items, such as three or more loci on a solid support.

[0083] As used herein, “distinct” refers to a unique value of molecular weight, mass shift or period that is different from every other value of molecular weight, mass shift or period in the measurement.

[0084] As used herein, “unambiguous” refers to the unique assignment of a particular oligonucleotide fragment according to the identity of its terminal base position and, in the event that a number of molecules are multiplexed, that the peak representing an oligonucleotide fragment can also be uniquely assigned to a particular molecule.

[0085] As used herein, the symbols M_(c), M_(T), M_(A) and M_(G) are average molecular weights in daltons of the nucleotides deoxycytidine, thymidine, deoxyadenosine and deoxyguanosine, respectively, or of analogs thereof. M_(avg), the average molecular weight of any given oligonucleotide is a function of the average molecular weights of each of the nucleotides comprising the oligonucleotide, the numbers c, t, a and g of each nucleotide present in the oligonucleotide, the length of the oligonucleotide n′ that is the sum of c, t, a and g, and the constant k that represents the mass of any other chemical groups on the molecule, such as terminal phosphates.

[0086] As used herein, N_(TOTAL) is the total number of possible base compositions for an oligonucleotide of length n′.

[0087] As used herein, “peak separation” or “minimum peak separation” S refers to the minimum value of the distance between consecutive peaks in a mass spectrum that resolves a large number of oligonucleotides having the same lengths but different molecular weights, i.e., different base compositions. The peak separation, which can be tailored by careful selection of the nucleotide analogs incorporated into the oligonucleotide and by a mass spectrometer of desired resolving power, is usually a positive integer greater than one, and typically a positive integer greater than or equal to 3. For two oligonucleotides having the same length n′ but different base compositions, their molecular weights will either correspond to the same peak if the molecular weights are identical, or to two peaks separated at least by a value equal to the peak separation.

[0088] As used herein, L is the maximum number of allowed oligonucleotide masses for a given nucleotide set. It is directly proportional to the oligonucleotide length n′ and the mass difference between the heaviest and lightest nucleotides in the set, and is inversely proportional to the peak separation.

[0089] As used herein, D refers to the average density of different base compositions per allowed mass value, given the set of all possible base compositions of an oligonucleotide of length n′.

[0090] As used herein, M_(heavy) refers to the mass of the heaviest nucleotide, nucleotide terminator or terminating base pair in daltons, depending on the specific embodiment of Forced Mass Modulation being described.

[0091] As used herein, M_(light) refers to the mass of the lightest nucleotide, nucleotide terminator or terminating base pair in daltons, depending on the specific embodiment of Forced Mass Modulation being described.

[0092] As used herein, M_(primer) is the mass of the primer in daltons.

[0093] As used herein, M_(obs)[n] is the observed mass of the sequencing reaction at the nth base position.

[0094] As used herein, M_(term)[n] refers to the mass in daltons of the nth terminating nucleotide.

[0095] As used herein, L′ is the theoretical upper limit on the number of sequences that be multiplexed in a single mass spectrum. L′ is directly proportional to the base periodicity P_(base), and is inversely proportional to the maximum mass shift S_(max).

[0096] As used herein, M_(duplex) is the mass in daltons of the fully duplex hairpin primer in the implementation of Forced Mass Modulation using pair-matched nucleotides.

[0097] As used herein, M_(ddM) is the mass in daltons of a dideoxy terminator that belongs to a set of mass-matched terminators.

[0098] As used herein, M_(targ)[n] is the mass of the nth nucleotide past the priming site in the 3′ to 5′ direction in the target sequence, i.e., the oligonucleotide whose sequence is being determined.

[0099] As used herein, “specifically hybridizes” refers to hybridization of a probe or primer only to a target sequence preferentially to a non-target sequence. Those of skill in the art are familiar with parameters that affect hybridization; such as temperature, probe or primer length and composition, buffer composition and salt concentration and can readily adjust these parameters to achieve specific hybridization of a nucleic acid to a target sequence.

[0100] As used herein, a biological sample refers to a sample of material obtained from or derived from biological material, such as, but are not limited to, body fluids, such blood, urine, cerebral spinal fluid and synovial fluid, tissues and organs. Derived from means that sample can be processed, such as by purification or isolation and/or amplification of nucleic acid molecules.

[0101] As used herein, a composition refers to any mixture. It may be a solution, a suspension, liquid, powder, a paste, aqueous, non-aqueous or any combination thereof.

[0102] As used herein, a combination refers to any association between two or among more items.

[0103] As used herein, “kit” refers to a package that contains a combination and optionally instructions and/or reagents and apparatus for use with the combination.

[0104] Forced Mass Modulation for Analysis of Nucleic Acid Molecules Time of Flight Analysis and Drawbacks Thereof.

[0105] While time-of-flight mass spectrometry offers a number of advantages over conventional techniques such as gel electrophoresis, the peculiar relationship between the masses of the bases in DNA complicates the analysis of complex mixtures of oligonucleotides by mass spectrometry. For a given oligonucleotide, the average molecular weight, M_(avg), is given by the following equation:

M _(avg) =k+cM _(c) +tM _(T) +aM _(A) +gM _(G)  (i)

[0106] where M_(c), M_(T), M_(A), M_(G) are the average molecular weights of each of the four nucleotide bases (cytosine, thymine, adenine, guanine) and c, t, a, g represent the number of each base present in the oligonucleotide. The term k is a constant representing the mass of any other chemical groups on the molecule, such as terminal phosphates. Rearranging equation (i) to give the average molecular weight as a function of the length of the oligonucleotide in bases yields

M _(avg) =k+n′M _(c) +t(M _(T) −M _(c))+a(M _(A) −M _(c))+g(M _(G−) M _(c))  (ii)

[0107] where n′, the oligonucleotide length, is defined as

n′=c+t+a+g

[0108] Substituting the masses of the naturally occurring bases in DNA (to one-tenth dalton):

M_(c)=289.2

M_(T)=304.2

M_(A)=313.2

M_(G)=329.2

[0109] into equation (ii) yields

M _(avg) =k+289.2n′+t(304.2−289.2)+a(313.2−289.2)+g(329.2−289.2),  (iii)

[0110] which can be simplified to

M _(avg) =k+289.2n′+15t+24a+40g  (iv)

[0111] Close inspection of equation (iv) reveals that it is almost always possible to find two oligonucleotides of the same length but of different base composition whose average masses differ by only one dalton. For example, all 7-mers having a base composition of A₂C₂ G₂ T have an average molecular weight of (2167.4+k), while all 7-mers with the base composition A₃CGT₂ have an average molecular weight of (2166.4+k). Since the following relation

(M _(c) +M _(G))=(M _(T) +M _(A))+1

[0112] is always true for the naturally occurring bases in DNA, simply replacing one C and one G in an oligonucleotide with one A and one T will produce a new oligonucleotide exactly one dalton lighter. Many other “single-dalton difference” relations, such as

4M _(A)=(M _(c) +M _(T)+2M _(G))+1

[0113] can readily be found for the naturally occurring bases.

[0114] Thus, the possibility always exists that two or more oligonucleotides of same length and different molecular weight (and, therefore, different base composition) will be too close in mass to be resolved by a time-of-flight instrument. Two oligonucleotides of same length but different molecular weight differ in base composition unless they are each composed of different nucleotide analogs, whereas two oligonucleotides of same length and same molecular weight can have either the same or different base compositions. This problem becomes increasingly severe with increasing oligonucleotide size, since the total number of possible base compositions, N_(TOTAL), scales as a cubic function of the oligonucleotide length n′, in bases: $\begin{matrix} {N_{TOTAL} = \frac{\left( {n^{\prime} + 1} \right)\left( {n^{\prime} + 2} \right)\left( {n^{\prime} + 3} \right)}{6}} & (v) \end{matrix}$

[0115] The use of time-of-flight mass spectrometry in sequencing applications also poses several potential problems. The great drawback of sequencing by the Sanger method is that the molecular weights of the Sanger reaction products can appear virtually anywhere on the mass axis depending on the particular sequence being examined. As a result, the absolute mass of any single Sanger fragment has to

[0116] be measured with sufficient accuracy to calculate its distance from the masses from the fragments above and below it. Thus, determination of the identity of a single base depends on the accuracy of two separate mass measurements. Any error in a determination mass of a single fragment affects the accuracy of two bases in the sequence.

[0117] For longer sequences (30-50 bases), it may not be possible to determine the mass difference between adjacent peaks with sufficient accuracy to unambiguously determine base identity. This is particularly a problem for the nucleotides A and T, which differ in mass only by nine daltons. The problem is addressed by resolving each of the four termination reactions in a separate mass spectrum. In this case each peak functions essentially as binary signal indicating the presence of a base at a particular position, much as in conventional electrophoretic sequencing. Using separate spectra, however, increases read accuracy but at the expense of increasing the number of required mass measurements by a factor of four.

[0118] It is possible to resolve two target sequences by the Sanger method in a single mass spectrum, provided that all products of the sequencing reactions have unique and resolvable masses, and multiplex methods using mass modified bases have been developed. But, where two or more reaction products have the same mass, then unambiguous reconstruction of the two target sequences is not possible (see, e.g., FIGS. 1c-e). In addition, there is no way to determine a priori which observed masses belong to a particular sequence. In practice, this means that multiplexed Sanger sequencing by mass spectrometry can be difficult. The methods provided herein resolve these problem and provide a way to determine a priori which masses are associated with extension of a particular primer.

[0119] Forced Mass Modulation

[0120] As noted above, Forced Mass Modulation refers to methods provided herein that permit unambiguous assignment of peak positions (or masses) to mass fragments of oligonucleotides according to their base position, base identity, and target sequence from which the fragments arose. The methods use deoxynucleotide analogs, modified nucleotide terminators, mass-labeled primers, mass-staggered primers and other such nucleotides, nucleic acids and analogs thereof to provide a means for deconvoluting complex mass spectra or output from other mass determining techniques. These methods permit deconvolution of highly multiplexed nucleic acid reaction mixtures for sequencing methods and detection methods that include a step of primer extension. In practicing these methods, primers are extended using mass-matched nucleotides and chain terminators (or in some embodiments mass where it is only necessary to detect incorporation (or the absence of incorporated) mass-matched terminators and optionally mass-matched chain extending nucleotides). Because the sequence and/or molecular mass of a primer is known, and the extended nucleotides have the same molecular mass, a periodicity in molecular mass that is a function of molecular weight of the selected mass matched nucleotide(s) results.

[0121] As described in more detail below, for sequencing reactions using chain terminators, the deviation from the periodicity results from incorporation of a chain terminator. The deviation is a function of the particular terminator incorporated. For detection methods, incorporation of a terminator will indicate the presence of a mutation (if the terminator is selected to pair with the first mutated nucleotide). Any shift from periodicity will indicate the presence of the mutation. These methods, thus provide a simple, reliable way to detect the presence of a mutations or target nucleotide(s) in a sequence and to sequence nucleic acids. Forced mass modulation can be used with any method, such as mass spectrometry and gel electrophoresis, that relies on molecular weight as an output. Mass spectrometry is exemplified herein.

[0122] The methods, designated Forced Mass Modulation methods, provided herein, are implemented by suitable selection of nucleotides and/or chain terminators, such as by the use of deoxynucleotide analogs, modified nucleotide terminators and mass-labeled primers in one or more reactions. Forced Mass Modulation can be used to simultaneously sequence or detect large numbers of oligonucleotides, such as, for example, twenty-five or more, with a high degree of resolution and accuracy. It can also be used to simplify the analysis of closely related sequence variants, as is required in the detection and scoring of nucleotide polymorphisms, including single nucleotide polymorphism (SNPs) and for other genotypical analyses. Forced Mass Modulation greatly improves the use of mass spectrometry for nucleic acid analyses. Nearly every application relies on mass measurements that can benefit in increased accuracy and in a reduction of the number of required spectra. Another advantage of Forced Mass Modulation is the number of different ways in which it can be implemented, allowing it to be tailored to particular experimental or instrumental limitations.

[0123] For example, compared to the conventional Sanger methods, Forced Mass Modulation, provides increased accuracy, simplified interpretation of mass data, and the ability to use a single mass spectrum for the unambiguous resolution of several distinct nucleic acid molecules. For mass spectrometry applications, the methods provide unambiguous assignment of peak positions of mass fragments of oligonucleotides according to their base position, base identity, and target sequence from which the fragments arose. Thus, the methods herein are advantageously used for multiplexing, in which a plurality of reactions are run in a single reaction (single pot). Forced Mass Modulation, exemplified with reference to sequencing methods, such as PROBE, can also be adapted to detection methods in which a primer is extended.

[0124] In Forced Mass Modulation in which a primer is extended with mass-matched nucleotides, for examples, the molecular weights of extended nucleic acid chains, such as sequencing reaction products, are constrained since all extension products from the same primer will have a molecular weight that differs either by the length of the extension and the chain terminator. As a result, the extension products assume a quasi-periodic distribution on the mass axis with a predetermined base periodicity. For sequencing, the base sequence itself is encoded in the pattern in which the observed mass distribution deviates from absolute regular periodicity. Since the base periodicity will always be known a priori, since the primer is known, each peak in the observed mass spectrum can be matched unambiguously to a unique nucleotide position in the target sequence. The initiating primers fix each set of nested fragments or extended products, and the use of mass-matched nucleotides creates the periodicity.

[0125] As demonstrated by the Examples below, the method is advantageous for numerous applications including sequencing and a variety of detection methods, including primer oligo base extension (PROBE) (see, e.g., U.S. Pat. No. 6,043,031; allowed U.S. application Ser. No. 09/287,679, now U.S. Pat. No. 6,258,538; and U.S. Pat. No. 6,235,478) that use mass spectrometry to distinguish between extended primers. If the base compositions of the target oligonucleotides are known a priori then it is possible to select a nucleotide set that produces oligonucleotide masses that are distinct and resolvable for any particular instrument or application.

[0126] Conversely, it is also possible to select a nucleotide set that restricts specific oligonucleotides to have the same mass, regardless of a change in base composition. The strategy of restricting specific oligonucleotides to have the same mass can be used to separate more than one oligonucleotide population of different lengths by restricting all oligonucleotides of a particular length to the same molecular weight, irrespective of differences in base composition.

[0127] The oligonucleotide analysis or sequencing in methods provided herein can be accomplished by one of several methods employed in the art for the synthesis, resolution and/or detection of nucleic acids. Depending on the embodiment implemented, modified nucleotides can be incorporated into the oligonucleotides by chemical (Oligonucleotides and Analogues: A Practical Approach, F. Eckstein, ed., IRL Press Oxford, 1991) or enzymatic (F. Sanger et al., Proc. Natl. Acad. Sci. USA 74:5463-67, 1977) synthesis. Extension products or truncated products of the oligonucleotides to be sequenced can be obtained using chemical (A. M. Maxam and W. Gilbert, Proc. Natl. Acad. Sci. USA 74:560-64, 1977) or enzymatic (F. Sanger et al., Proc. Natl. Acad. Sci. USA 74:5463-67, 1977) methods.

[0128] For the resolution and detection of target nucleic acids any mass determination method, such as, but are not limited to, chromatography, gel electrophoresis, capillary zone electrophoresis and mass spectrometry, is used. Mass spectrometric formats, include, but are not limited to, are matrix assisted laser desorption ionization (MALDI), electrospray (ES), ion cyclotron resonance (ICR) and Fourier Transform. For ES, the samples, dissolved in water or in a volatile buffer, are injected either continuously or discontinuously into an atmospheric pressure ionization interface (API) and then mass analyzed by a quadrupole. The generation of multiple ion peaks which can be obtained using ES mass spectrometry can increase the accuracy of the mass determination. Even more detailed information on the specific structure can be obtained using an MS/MS quadrupole configuration

[0129] In MALDI mass spectrometry, various mass analyzers can be used, e.g., magnetic sector/magnetic deflection instruments in single or triple quadrupole mode (MS/MS), Fourier transform and time-of-flight (TOF) configurations as is known in the art of mass spectrometry. For the desorption/ionization process, numerous matrix/laser combinations can be used. Ion-trap and reflectron configurations can also be employed.

[0130] Pair-matched Nucleotide-based Methods

[0131] Forced Mass Modulation can be implemented using a deoxynucleotide set in which the mass of each base pair is identical, termed a pair-matched nucleotide set. A pair-matched nucleotide set can easily be formed, for example, by replacing dG (329.2 Da) in the set of naturally occurring nucleotides with 7-deaza-dG (328.2 Da). This forces the mass of each base pair to be 617.4 daltons:

(dA+dT)=(313.2+304.2)=617.4 Da

(dC+7-deaza-dG)=(289.2+328.2)=617.4 Da

[0132] Many other pair-matched sets are possible using available nucleotide analogs. For this embodiment, the target DNA sequence can be composed entirely of the pair-matched nucleotide set. This can be accomplished by amplifying the target DNA sequence by PCR using the pair-matched nucleotide set prior to the sequencing reaction.

[0133] A further requirement for this embodiment of Forced Mass Modulation is that the mass of each terminating base pair is unique and resolvable. The standard dideoxy terminators therefore cannot be used with the pair-matched nucleotide set described above, because the masses of all terminating base pairs are identical at 601.4 daltons, except ddG:dC, which is 602.4 daltons. For the sake of clarity in this example, it is assumed that a set of mass-matched terminators is available (“ddM,” defined as set of chain-terminating nucleotides that have exactly the same molecular weight ddA=ddC=ddG=ddT). If the mass of ddM is arbitrarily chosen to be 500 daltons, then the masses of the terminating base pairs are as follows: Terminating Base Pair Mass (Da) ddM: dC 789.2 ddM: dT 804.2 ddM: dA 813.2 ddM: 7-deaza-dG 828.2

[0134] In practice it is also possible to implement Forced Mass Modulation using a set of terminators that have different masses, this is discussed in detail below.

[0135] Exemplary embodiments in which the mass shift is obtained using pair-matched nucleotides, where the mass of each nucleotide base-pair is selected so that the masses of all pairs are identical, are described in Example 4. In one embodiment thereof, the following steps are performed: (i) the target nucleic acid is copied or amplified by a method such as PCR in the presence of the pair-matched nucleotide set prior to the sequencing or detection reaction; (ii) the target nucleic acid is denatured, and a partially duplex hairpin primer is annealed and ligated to the single-stranded template; (iii) the primer is extended in the presence of chain terminating nucleotides and pair-matched nucleotides to produce extension products; (iv) the masses of the extension products follow a periodic distribution that is determined by the mass of the pair-matched nucleotide set, and, (v) the target nucleic acid is detected by virtue of its molecular weight or its sequence is determined from the mass shift of each extension product from its corresponding periodic reference mass.

[0136] In embodiments described above, the extending bases are pair matched. the mass of each terminating base pair is unique and resolvable, so that the mass shifts corresponding to each terminating base pair are unique. The nucleotide terminators are optionally mass-matched or can be of distinct masses as long as distinct values of mass shift are obtained for each terminating base pair.

[0137] In another embodiment, the extension products are treated to produce blunt-ended double-stranded extension products by methods known to those of skill in the art, such as the use of single-strand specific nucleases. In an aspect of this embodiment, a plurality of target nucleic acids can be multiplexed in a single reaction by annealing each target nucleic acid to a primer of distinct molecular weight. The primers can be selected to differ in molecular weight by a value that is greater than the maximum mass shift, i.e., the difference in molecular weight between the heaviest and the lightest nucleotide terminating base pairs. Since double stranded nucleic acid can be analyzed, the effective sequence read is halved relative to the embodiment employing mass-matched nucleotides, but the number of molecules that can be multiplexed is doubled, due to the increase in period (the value of the mass of a base pair, rather than a single mass-matched nucleotide). In exemplary embodiments, about 14 to about 50 sequences are multiplexed. In detection embodiments, about 50 or more molecules can be simultaneously multiplexed since only a single terminating base pair is added in the extension reaction.

[0138] In another embodiment, the chain termination reactions can each be carried out separately using a standard nucleotide terminator, pair-matched nucleotides, and mass-labeled primers, if modified nucleotide terminators which are either mass-matched or provide distinct values of mass shift for each terminating base pair are not available. The reactions can be pooled prior to detection or sequence analysis. In one embodiment, the mass-labeled primers can have distinct values of molecular weight that give rise to unique values of mass shift or positional mass difference for each terminating base.

[0139] Optimizing the Mass Spectrometric Analysis of Oligonucleotide Mixtures

[0140] In another method provided herein, nucleotide analogs are used to restrict the possible values of molecular weights that an oligonucleotide can possess relative to other oligonucleotides of the same length. The nucleotide analogs can be incorporated into the oligonucleotides using any suitable method, such as automated DNA synthesis (Oligonucleotides and Analogues: A Practical Approach, F. Eckstein, ed., IRL Press Oxford, 1991) or by enzymatic replication using a polymerase and the requisite nucleotides and nucleotide analogs.

[0141] For example, any two oligonucleotides with the same length n′ with different base compositions can either 1) have exactly the same average molecular weight, or 2) have molecular weights no closer than a minimum value called the peak separation. In most cases, the peak separation will be a positive integer greater than one, but fractional values are theoretically possible.

[0142] To illustrate an exemplary implementation of this method, the average molecular weight of the nucleotide analog 7-deaza-dG (328.2 daltons) can be substituted for MG, into equation (ii) above, which defines M_(avg) as a function of the length “n” of the oligonucleotides in bases, as follows:

M _(avg) =k+n′M _(c) +t(M_(T) −M _(c))+a(M _(A) −M _(c))+g(M _(G) −M _(c)), where

[0143] M_(c), M_(T), M_(A), M_(G) are the average molecular weights of each of the four nucleotide bases (cytosine, thymine, adenine, guanine); c, t, a, g represent the number of each base present in the oligonucleotide, the sum thereof, i.e., c+t+a+g=n′, the total oligonucleotide length in bases; and the term k is a constant representing the mass of any other chemical groups on the molecule, such as terminal phosphates. Substituting the masses of the naturally occurring bases dC, dT and dA in DNA (to one-tenth dalton), and of 7-deaza-dG,

M_(c)=289.2

M_(T)=304.2

M_(A)=313.2

M_(G)=328.2

[0144] and following simplification, the equation reduces to:

M _(avg) =k+289.2n′+15t+24a+39g

[0145] Extracting the common factor from the last three terms yields

M _(avg) =k+289.2n′+(5t+8a+13g)×3  (iv)

[0146] In this example, the minimum peak separation is three daltons. It is not possible to identify or detect two oligonucleotides of the same length with different molecular weights that are closer than three daltons. Oligonucleotides with average masses closer than three daltons the oligonucleotides are detected if they are of different lengths.

[0147] As a second example, MT can be substituted with the molecular weight of a hypothetical nucleotide analog whose mass is 305.2 into equation (ii), yielding

M _(avg) =k+289.2n′+16t+24a+40g

[0148] Extracting the common integer factor from the last three terms yields

M _(avg) =k+289.2n′+(2t+3a+5g)×8  (vii)

[0149] for a minimum peak separation of eight daltons. Thus, appropriate selection of nucleotide analogs permits construction of nucleotide sets that provides sufficient peak separation for adequate resolution by mass, such as in a time-of-flight mass spectrometer. The trade-off for a greater peak separation is a greater number of base compositions that have exactly the same mass for a given oligonucleotide length. The maximum number of allowed oligonucleotide masses, L, for a given nucleotide set, is given by $\begin{matrix} {{L = {\frac{n^{\prime}\left( {M_{heavy} - M_{light}} \right)}{S} + 1}},} & ({viii}) \end{matrix}$

[0150] where n′ is the oligonucleotide length in bases, S is the peak separation, M_(light) the mass of the lightest nucleotide in the set, M_(heavy) is the mass of the heaviest nucleotide in the set. The number of allowed oligonucleotide masses scales in direct proportion to the base length and inversely with the peak separation, but not all possible mass values will be represented for a given oligonucleotide length, particularly for small n. The average density of different base compositions per allowed mass value, D, can be obtained by dividing equation (v) by (viii) ${D = \frac{N_{TOTAL}}{L}},$

[0151] which expands into $\begin{matrix} {D = \frac{{S\left( {n^{\prime} + 1} \right)}\left( {n^{\prime} + {2\left( {n^{\prime} + 3} \right)}} \right.}{6\left( {{n^{\prime}\left( {M_{G} - M_{C}} \right)} + S} \right)}} & ({ix}) \end{matrix}$

[0152] using a typical nucleotide set with G as the heaviest base and C as the lightest. The density function scales in direct proportion to the peak separation and as a quadratic function of the oligonucleotide length in bases. In practice, the average density of base compositions per allowed mass value predicated by equation (ix) will be somewhat lower than the actual density of base compositions per observed mass value, because not all allowed masses will always be represented. The Examples describe implementation of the methods for sequencing.

[0153] System and Software method for Force Mass Modulation

[0154] Also provided are systems that automate the methods for determining a nucleotide sequence of a target nucleic acid or the detection methods provided herein using a computer programmed for identifying the sequence or target nucleic acid identity based upon the methods provided herein. The methods herein can be implemented, for example, by use of the following computer systems and using the following calculations, systems and methods.

[0155] An exemplary automated testing system contains a nucleic acid workstation that includes an analytical instrument, such as a gel electrophoresis apparatus or a mass spectrometer or other instrument for determining the mass of a nucleic acid molecule in a sample, and a computer capable of communicating with the analytical instrument (see, e.g., copending U.S. application Ser. Nos. 09/285,481, 09/663,968 and 09/836,629; see, also International PCT application No. WO 00/60361 for exemplary automated systems). In an exemplary embodiment the computer is an IBM compatible computer system that communicates with the instrument using a known communication standard such as a parallel or serial interface.

[0156] For example, systems for analysis of nucleic acid samples are provided. The systems include a processing stations that performs a forced mass modulation chain extension reaction; a robotic system that transports the resulting products from the processing station to a mass measuring station, where the masses of the products of the reaction are determined; and a data analysis system, such as a computer programmed to identify nucleotides using forced mass modulation data, that processes the data from the mass measuring station to identify a nucleotide or plurality thereof in a sample or plurality thereof. The system can also include a control system that determines when processing at each station is complete and, in response, moves the sample to the next test station, and continuously processes samples one after another until the control system receives a stop instruction.

[0157] The computer can be part of the instrument or another system component or it can be at a remote location. A computer system located at a site distant from the instrument can communicate with the instrument, for example, through a wide area network or local area communication network or other suitable communication network. The system with the computer is programmed to automatically carry out steps of the methods herein and the requisite calculations. For embodiments that use mass-matched deoxyriboucleotides, a user enters the primer sequence or primer mass, the periodic reference mass and mass of an individual mass-matched deoxyonucleotide. These data can be directly entered by the user from a keyboard or from other computers or computer systems linked by network connection, or on removable storage medium such as a CD-ROM, minidisk (MD), DVD, floppy disk or other suitable storage medium. Next the user causes execution software that operates the system in which the mass spectrum of the extension products is generated. The Forced Mass Modulation software performs the steps of obtaining the masses of the fragments generated by the sequencing reaction and measured by the analytical instrument, and determining the identity of a nucleotide at any base position or the positional mass difference. The identity of the nucleotide at each base position is determined by comparing the calculated M_(diff)[n] values to a database of previously calculated values of M_(diff) for each of the chain terminating nucleotides.

M _(diff) [n]=M _(obs) [n]−M _(PR) [n],

[0158] where:

M _(PR) [n]=(M _(primer) +M _(light))+(n−1) P _(base,)  (i)

[0159] in which n is the base position, M_(PR)[N] is the n^(th) periodic reference mass, M_(primer) is the mass of the primer, M_(light) is the mass of the lightest nucleotide terminator and P_(base) is the base periodicity in daltons. The observed masses of the sequencing reaction products are given by the following equation:

M _(obs) [n]=M _(primer)+(n−1) P _(base) +M _(term) [n],  (ii)

[0160] where n is the base position, M_(obs)[n] is the nth observed mass, P_(base) is the base periodicity, and M_(term)[n] is the mass of the nth terminating nucleotide in daltons. The positional mass differences for the sequence can be obtained by subtracting equation (i) from equation (ii) and evaluating at every base position n:

[0161] where M_(diff)[n] is the n^(th) positional mass difference. This relation simplifies to:

M _(diff) [n]=M _(term) [n]−M _(light).  (iii)

[0162] Hence, the periodicity is determined by the mass of the mass-matched nucleotide and the shift is the difference in location of a peak resulting from the chain terminator. For example, in FIG. 2, the lightest terminator is ddC, and the differential is 0 for C, 40 for G, 34 for A, 15 for T. The selected mass matched nucleotide has a mass of 310 Da. The primer in FIG. 2a has a mass of 3327 Da and the first peak would be at 3600 if the first nucleotide in the extension product were C (0 shift). Since the first peak is at 3640, the shift is 40 Da. Therefore the first nucleotide is G, corresponding to a shift from the periodicity of 310 Da generated by the mass-matched nucleotides.

[0163] Detection Methods

[0164] The methods herein may be used with any method for detection of nucleic acids based on molecular mass known to those of skill in the art, particularly methods in which a primer is extended. Such methods are modified by extending using mass matched nucleotides and/or chain terminators in extension reactions. Alternatively, or additionally, amplification reactions may be performed using mass-matched nucleotides or pair-matched sets of nucleotides. These methods can be readily multiplexed using the methods and nucleic acid molecules provided herein.

[0165] Detection methods and protocols, including those that rely on mass spectrometry (see, e.g., U.S. Pat. Nos. 6,194,144; 6,225,450; 5,691,141; 5,547,835; 6,238,871; 5,605,798; 6,043,031; 6,197,498; 6,235,478; 6,221,601; 6,221,605; International PCT application No. WO 99/31273, International PCT application No. WO 98/20019), can be modified for use with the methods herein by using mass-matched nucleotides for extension or pair matched duplexes for hybridization reactions.

[0166] Among the methods of analysis herein are those involving the primer oligo base extension (PROBE) reaction with mass spectrometry for detection. In such reactions, the primer will be extended by mass-matched nucleotides. The methods herein are designed for multiplexing so that a plurality of different primers can be extended at different loci in the same reaction. The PROBE method uses a single detection primer followed by an oligonucleotide extension step to give products, which can be readily resolved by mass spectrometry, and, in particular, MALDI-TOF mass spectrometry. The products differ in length depending on the presence or absence of a polymorphism. In this method, a detection primer anneals adjacent to the site of a variable nucleotide or sequence of nucleotides and the primer is extended using a DNA polymerase in the presence of one or more dideoxy NTPs and, optionally, one or more deoxy NTPs. The resulting products are resolved by MALDI-TOF mass spectrometry. The mass of the products as measured by MALDI-TOF mass spectrometry makes possible the determination of the nucleotide(s) present at the variable site. Use of primers containing mass-matched bases increases the resolving power of the reaction and permit simultaneous detection of a plurality of mutations (polymorphisms).

[0167] These methods can be automated (see, e.g., copending U.S. application Ser. No. 09/285,481 and published International PCT application No. PCT/US00/08111, which describes an automated process line) and performed in a system that includes a computer programmed for analysis of the mass data as described above.

[0168] The analyses can be performed on chip based formats in which the target nucleic acids or primers are linked to a solid support, such as a silicon or silicon-coated substrate, preferably in the form of an array. Generally, when analyses are performed using mass spectrometry, particularly MALDI, small nanoliter volumes of sample are loaded on, such that the resulting spot is about, or smaller than, the size of the laser spot. It has been found that when this is achieved, the results from the mass spectrometric analysis are quantitative. The area under the signals in the resulting mass spectra are proportional to concentration (when normalized and corrected for background). Methods for preparing and using such chips are described in U.S. Pat. No. 6,024,925, co-pending U.S. application Ser. Nos. 08/786,988, 09/364,774, 09/371,150 and 09/297,575; see, also U.S. application Serial No. PCT/US97/20195, which published as WO 98/20020. Chips and kits for performing these analyses are commercially available from SEQUENOM under the trademark MassARRAY. MassArray relies on the fidelity of the enzymatic primer extension reactions combined with the miniaturized array and MALDI-TOF (Matrix-Assisted Laser Desorption Ionization-Time of Flight) mass spectrometry to deliver results rapidly. It accurately distinguishes single base changes in the size of DNA fragments associated with genetic variants without tags.

[0169] The following Examples are included for illustrative purposes only and are not intended to limit the scope of the invention.

EXAMPLE 1

[0170] Forced Mass Modulation using Mass-Matched Deoxynucleotides

[0171] For this implementation, a set of nucleotide analogs for the four bases in DNA are selected (Amersham Pharmacia Biotech) such that each base has exactly the same molecular weight, termed a mass-matched deoxynucleotide set. This is achieved by judiciously choosing chemical modifiers of the existing bases or by the using a universal base analog such as deoxyinosine, which can form base pairs with more than one of the natural bases. For this example, the mass of each deoxynucleotide (“dN”) in the mass-matched set has the arbitrarily selected value of 310 daltons, but any other value suffices. The sequencing reaction is performed as follows: 1) a primer is annealed to the target to be sequenced; 2) the resulting structure is subjected to an extension reaction using a suitable polymerase in the presence of the mass-matched nucleotide set and the four standard dideoxynucleotide terminators. The products and molecular masses of such a reaction are shown with a simulated mass spectrum in FIG. 2a. The base periodicity is the mass of dN, or 310 daltons. The identity of a nucleotide at any base position is given by the positional mass difference, defined as the distance in daltons between the observed peak and the nearest periodic reference mass, which occurs every 310 daltons. In this example, the first periodic reference mass is defined as the (primer mass+ddC), or (3327+273)=3600 daltons. The second periodic reference mass would be 3600 plus the base periodicity or (3600+310)=3910, and so on. Expressed in terms of the base position n:

M _(PR) [n]=(M _(primer) +M _(light))+(n−1) P _(base),  (i)

[0172] where n is the base position, M_(PR) [n] is the n ^(th) periodic reference mass, M_(primer) is the mass of the primer, M_(light) is the mass of the lightest nucleotide terminator and P_(base) is the base periodicity in daltons. The observed masses of the sequencing reaction products are given by the following equation:

M _(obs) [n]=M _(primer)+(n−1) P _(base) +M _(term) [n],  (ii)

[0173] where n is the base position, M_(obs)[n] is the n^(th) observed mass, P_(base) is the base periodicity, and M_(term)[n] is the mass of the n^(th) terminating nucleotide in daltons. The positional mass differences for the sequence can be obtained by subtracting equation (i) from equation (ii) and evaluating at every base position n:

M _(diff) [n]=M _(obs) [n]−M _(PR) [n],

[0174] where M_(diff)[n] is the n^(th) positional mass difference. This relation simplifies to:

M _(diff) [n]=M _(term) [n]−M _(light).  (iii)

[0175] Inspection of equation (iii) reveals that M_(diff) can only take on four distinct values, each corresponding to a different nucleotide terminator:

M _(diff) [“ddC”]=(273.2−273.2)=0

M _(diff) [“ddT”]=(288.2−273.2)=15

M _(diff) [“ddA”]=(297.2−273.2)=24

M _(diff) [“ddG”]=(313.2−273.2)=40.

[0176] Hence, the identity of the nucleotide at every base position in the target sequence can be determined by comparing each calculated positional mass difference with the values in the table above. Since the values that M_(diff) can assume depend only on the choice of nucleotide terminators used in the sequencing reaction, it is possible to tailor the positional mass differences so that they are resolvable for any particular mass spectrometer. For example, replacing the terminator ddT with its analog 5-bromo-dideoxyuridine (353.1 daltons) yields a positional mass difference of (353.1−273.2)=79.9 Da for termination at T positions in the target sequence. This type of nucleotide substitution can be particularly valuable for lower-resolution mass spectrometers, as it possible to maintain the sequence read accuracy without requiring any additional mass spectra.

[0177] Further inspection of equation (iii) reveals that each observed mass value can be at most 40 daltons heavier than the nearest periodic reference mass. This limit is termed the maximum mass shift and is defined as the mass difference between the heaviest nucleotide terminator and the lightest. Resolving a second target sequence by Forced Mass Modulation with the standard dideoxy terminators is possible in a single spectrum so long as the primer for the second sequence is at least 40 daltons heavier (the maximum mass shift) than the primer for the first sequence, thus insuring that the peaks for each sequence never overlap in mass.

[0178] In practice, it is recommended most mass spectrometric formats that the second primer is at least about 60 daltons heavier than the first primer, as each observed peak will have a finite width. FIG. 2b shows a target second sequence resolved on the same mass spectrum shown in FIG. 2a, using a primer heavier by 77 daltons. The peaks corresponding to the reaction products from the first target sequence can fall within the shaded regions of the spectrum in FIG. 2b, which can never intersect peaks from the second target sequence. Unambiguous resolution of both sequences is possible in this arrangement because each peak can be uniquely assigned to a nucleotide, a base position, and a target sequence. This method is designated Mass Spectrum Division Multiplexing herein, and it is implemented using mass-staggered primers. FIG. 3 shows four different sequences resolved in a single spectrum using a set of mass-staggered primers that are separated in mass by integer multiples of 77 daltons (77, 154, and 231 daltons).

[0179] The theoretical upper limit on the number of sequences that can be multiplexed in a single mass spectrum is given by the following equation: $\begin{matrix} {L^{\prime} = \frac{P_{base}}{S_{\max},}} & ({iv}) \end{matrix}$

[0180] where L′ is the upper limit, P_(base) is the base periodicity, and S_(max) is the maximum mass shift in daltons. For the nucleotide set and terminators used in this example, L=(310/40)=7.75, or approximately seven. Increasing the number of sequences that can be multiplexed in a single spectrum, can be achieved by implementing one or both of an increase in the base periodicity, and a reduction of the maximum mass shift. The base periodicity can be increased by choosing a mass-matched nucleotide set that has a higher molecular weight for dN. It is simpler to lower the maximum mass shift by careful use of the nucleotide terminators and their analogs. For example, if the sequencing reactions were performed using only the terminators ddC, ddT, and ddA, then the maximum mass shift becomes (mass of ddA−mass of ddC)=(297-273)=24 Da. In this case the upper limit on the number of sequences that can be multiplexed is L=(310/24)=12.92, or approximately twelve. In situations where complete sequence information is not required, such as diagnostic sequencing, a great reduction in the number of required spectra can be realized by using fewer than four nucleotide terminators. If the sequencing reaction is performed using only a single nucleotide terminator, the maximum mass shift becomes identically zero, and the number of sequences that can be multiplexed in a single spectrum is limited only by the absolute resolution of the mass spectrometer in question. If a given mass spectrometer has an absolute resolution of 12 Da in the mass range of the sequencing reaction products, then the maximum number of sequences that can be multiplexed is given by L=(310/12)=25.83, or approximately twenty-five.

EXAMPLE 2

[0181] Forced Mass Modulation using Pair-Matched Deoxynucleotides

[0182] Implementation of Forced Mass Modulation using pair-matched nucleotides is shown in FIG. 4. The basic requirement for this method is that the sequencing reaction products can be analyzed as double-stranded structures. Briefly, the steps in the reaction are as follows: 1) A partially duplex hairpin primer with a 3′ overhang and a 5′ phosphate group is annealed and ligated to the single stranded target sequence. 2) The resulting partially duplex structure is subjected to a sequencing reaction using the pair-matched nucleotide set described above along with the set of mass-matched terminators (ddM). 3) The products from the sequencing reaction are exposed to a strict single strand-specific nuclease that results in the production of blunt-ended hairpin structures ready for analysis by mass spectrometry. FIG. 5 shows the products and molecular masses of the nuclease digestion along with a simulated mass spectrum.

[0183] Because the reaction products are double-stranded, they are forced to assume a quasi-periodic distribution with a base periodicity of 617.4 daltons. The shaded regions on the spectrum shown in FIG. 5 indicate the allowed mass ranges that can be occupied by the reaction products. The first periodic reference mass is at 10360 Da, which is the mass of the fully duplex hairpin primer plus a ddM:dC base pair. Expressing the periodic reference masses in terms of the base position n yields:

M _(PR) [n]=(M _(duplex) +M _(light) +Mdd _(M))+(n−1)×P _(base)  (x)

[0184] Where M_(PR)[n] is the n^(th) periodic reference mass, M_(duplex) is the mass of the fully duplex primer, M_(light) is the mass of the lightest deoxynucleotide in the target, P_(base) is the base periodicity, and Mdd_(M) is the mass of ddM in daltons. The observed masses of the sequencing reaction products are given by the following equation:

M _(obs) [n]=M _(duplex) +Mdd _(M)+(n−1)×P _(base) +M _(targ) [n],  (xi)

[0185] where n is the base position, M_(ob)[n] is the nth observed mass, and M_(targ)[n] is the mass of the n^(th) nucleotide in the target sequence past the priming site in the 3′->5′ direction.

[0186] In contrast to the mass-matched nucleotide set implementation that provides the sequence complementary to the template strand read in the 5′->3′ direction, the pair-matched nucleotide set implementation described herein directly reads the template strand in the 3′->5′ direction. The positional mass differences for this implementation are the same as those in Example 1, except that the mass difference corresponding to a termination on dG is 39 as opposed to 40 daltons, because 7-deaza-dG is exactly one dalton lighter than dG. Since double stranded DNA can be analyzed for this method to work, the effective sequence read length is halved, although the number of sequences that can be multiplexed is doubled, due to the increase in the base periodicity.

[0187] As a demonstration of Forced Mass Modulation implemented without using mass-matched terminators, the positional mass differences for the above example using the following set of nucleotide terminators is calculated as follows: Termi- Base Mass of nator Nucleotide Analog Mass Pairing Base Pair T 5-Bromo-dideoxyuridine 353.1 5-Br-ddU: dA 666.3 C 5-Methyl-dideoxycytidine 287.2 5-Me-ddC: 615.4 7-deaza-dG A Dideoxyadenosine 297.2 ddA: dT 601.4 G Dideoxyinosine 298.2 ddI: dC 587.4

[0188] The positional mass difference at every base position is given by:

M _(diff) [n]=M _(pair) [n]−M _(lightest),  (xii)

[0189] where M_(diff)[n] is the nth positional mass difference, M_(pair)[n] is the mass of n^(th) terminating base pair, and M_(lightest) is the mass of the lightest terminating base pair in daltons. Substituting in the values from the table above yields:

M _(diff) [“G”]=(587.4−587-4)=0

M_(diff) [“A”]=(601.4−587-4) 14

M _(diff) [“C”]=(615.4−587.4)=28

M _(diff) [“T”]=(666.3−587-4)=78.9

[0190] Since each terminating base pair has a unique positional mass difference, the base sequence can be determined unambiguously. The maximum mass shift in this case is 78.9 daltons. When choosing a set of terminating nucleotides it is important to select the set such that the positional mass difference for each base termination is distinct and resolvable by mass.

[0191] If modified nucleotide terminators are not used, it is still possible to implement Forced Mass Modulation by carrying out each of the four termination reactions separately using mass-labeled primers rather than modified terminators, combining all reaction products, and then obtaining a mass spectrum. In order to produce the same positional mass differences as shown in Example 1, using a set of pair-matched nucleotides and the standard dideoxy terminators, the following primer mass shifts are required: Termination Reaction Primer Mass C “reference” primer T reference primer + 15 Da A reference primer + 24 Da G reference primer + 39 Da

[0192] This method is essentially equivalent to multiplexing four single-nucleotide sequencing reactions in the same spectrum, except that all the sequencing products originate from the same priming site but terminate on different nucleotides.

EXAMPLE 3

[0193] Forced Mass Modulation in the Detection and Scoring of Single Nucleotide Polymorphisms

[0194] Forced Mass Modulation can be used to simplify the analysis of closely related sequence variants, as is required in the detection and scoring of single nucleotide polymorphisms. FIG. 6 shows three sequence variants that differ from each other only at a single base position sequenced by a conventional Sanger reaction. The mass distribution of the reaction products is so complex that it can be uninterpretable, even if the base sequences of the variants are known a priori.

[0195]FIG. 7 shows the same three variants sequenced by Forced Mass Modulation using mass-matched deoxynucleotides (dN=310 Da) and the standard dideoxy terminators. The positions and identities of the single-nucleotide changes are immediately apparent from the mass spectrum. Since the masses of the sequencing reaction products are constrained to fall within the shaded regions of the spectrum in FIG. 7b, it is possible to multiplex other sequences on the same spectrum.

EXAMPLE 4

[0196] Base Composition Density Distributions for the Total Set of possible 7-base Oligonucleotides

[0197] For this implementation, three sets of 7-base oligonucleotides comprising all possible base compositions for a 7-base oligonucleotide can be obtained; the first set comprising the four natural bases (dA, dG, dC and dT), the second set comprising three of the natural bases (dA, dC and dT) and the nucleotide analog 7-deaza-deoxyguanosine (7-deaza-dG) substituted for dG, and the third set comprising three of the natural bases (dA, dG and dC) and the nucleotide analog deutero-deoxythymine (deutero-dT) substituted for dT. FIG. 8 shows the actual base composition density distributions for the total set of possible 7-base oligonucleotides using the three different nucleotide sets. Note that for the set of naturally occurring bases (FIG. 8a), nearly every base composition has its own distinct mass value, but most of these mass values are spaced only one dalton from each other. Increasing the peak separation to three daltons by substitution of dG with 7-deaza-dG (FIG. 8b) markedly increases the average number of base compositions per observed mass, particularly for those masses in the center of the range, but any two oligonucleotides of the same length with different molecular weights will have to be separated by at least three daltons. Similarly, substitution of dT with deutero-dT (FIG. 8c) gives a minimum peak separation between oligonucleotides having the same length but different molecular weights of eight daltons. The trade-off for a greater peak separation is a greater number of oligonucleotides that have exactly the same mass for a given oligonucleotide length.

[0198] Since modifications will be apparent to those of skill in this art, it is intended that this invention be limited only by the scope of the appended claims.

1 1 1 18 DNA Artificial Sequence Oligonucleotide target sequence 1 tgaatgtaga cctcacgt 18 

What is claimed is:
 1. A method for identifying nucleotides at one or more base positions in a plurality of target nucleic acids molecules, comprising: synthesizing extension products of the target nucleic acid in the presence of chain terminating nucleotides and mass-matched nucleotides; determining the mass of each extension product; and calculating a mass shift from a period for the mass of each extension product, whereby the nucleotides in the target nucleic acid molecules are identified by determining the nucleotide that corresponds to each mass shift.
 2. The method of claim 1, wherein the mass-matched deoxy-nucleotides are identical.
 3. The method of claim 1, wherein a mass-matched deoxy-nucleotide is deoxyinosine, 5-nitroindole, 3-nitropyrrole, 3-methyl 7-propynyl isocarbostyril, 5-methyl isocarbostyril or 3-methyl isocarbostyril.
 4. A kit for determining the sequence of a target nucleic acid, comprising mass-matched nucleotides.
 5. A kit for determining the sequence of a target nucleic acid, comprising pair-matched nucleotides and mass-matched chain terminating nucleotides.
 6. A kit for determining the sequence of a target nucleic acid, comprising pair-matched nucleotides and chain terminating nucleotides that form base pairs of distinct molecular weight, and optionally including instructions for sequencing using these reagents.
 7. A kit for determining the sequence of a target nucleic acid, comprising pair-matched nucleotides and mass-labeled primers, and optionally including instructions for sequencing using these reagents.
 8. A computer-based method for identifying nucleotide or nucleotides at one or more base positions in a target nucleic acid molecule or plurality thereof, comprising: a) entering a computer, a primer sequence or primer mass, a mass of an individual mass-matched deoxyonucleotide and the identify of chain terminators used; b) entering the masses of the fragments generated by a primer extension reaction, wherein the primer is extended by mass-matched deoxynucleotides; c) determining P_(base), wherein P_(base) is the base periodicity in daltons; d) calculating M_(diff)[n] for each nucleotide base to be identified, wherein: M _(diff) [n]=M _(obs) [n]−M _(PR) [n]; M _(PR) [n]=(M _(primer) +M _(light))+(n−1) P _(base);  M_(obs)[n] is the observed peak;  where: n is the base position; M_(PR)[n] is the nth periodic reference mass; M_(primer) is the mass of the primer; M_(light) is the mass of the lightest nucleotide terminator; and e) determining the identity of a nucleotide at any base position or the positional mass difference by determining M_(diff)[n] and comparing it to a database of previously calculated values of M_(diff) for each of the chain terminating nucleotides.
 9. A system for high throughput analysis of nucleic acid samples, comprising: a processing station that performs a chain extension reaction, in the presence of mass-matched nucleic nucleotides, on a nucleic acid sample in a reaction mixture; a robotic system that transports the resulting products from the processing station to a mass measuring station, wherein the masses of the products of the reaction are determined; and a data analysis system that processes the data from the mass measuring station and that is programmed to perform the method of claim 46 to identify a nucleotide or nucleotides at one or more base positions in nucleic acid molecule in the sample.
 10. The system of claim 10, further comprising a control system that determines when processing at each station is complete and, in response, moves the sample to the next test station, and continuously processes samples one after another until the control system receives a stop instruction.
 11. The system of claim 10, wherein the mass measuring station is a mass spectrometer. 